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Introduction 




Internal (deep dielectric) charging 

High energy (>100 keV) electrons penetrate 

spacecraft walls and accumulate in dielectrics 
or isolated conductors 

Threat environment is energetic electrons with 
sufficient flux to charge circuit boards, cable 
insulation, and ungrounded metal faster than 
charge can dissipate 

Accumulating charge density generates electric 
fields in excess of breakdown strength 
resulting in electrostatic discharge 

System impact is material damage, discharge 
currents inside of spacecraft Faraday cage on 
or near critical circuitry, and RF noise 

Overview 

• Internal charging physics, model formulation 

• Examples of two engineering models 
ESA DICTAT 

NASANUMIT 

• Other internal charging models 


PMMA (acrylic) charged by ~2 to 5 MeV electrons 
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Day in solar rotation period 
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GOES Solar Cycle 21 Internal Charging Anomalies (GEO) 



smoothed sunspot number 
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• Black: 


GOES phantom commands 


[Wrenn et al. 2002 ] 


• Red: 

• Amber: 

• Green: 

• White: 


2-day fluence (F2) > 2 MeV electrons 

F2 > 10 9 e/cm 2 -sr 
10 9 > F2 > 10 8 e/cm 2 -sr 
F2 < 10 8 e /cm 2 -sr 
no data 
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Internal Charging: Physics 


0.5 < a < 1.0 


Electron 

Beam 



Conductor CD Ampere meter. Arrow shows 
^ positive current direction. 


[Jun et al. 2007] 


Solution to Poisson, continuity equation involves two problems: 

• Radiation (electron) penetration with charge and energy deposition in material 

• Electrostatic solution of fields from charge distribution in insulator 
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Internal Charging: 


Computation of radiation dose and 
charge deposition as function of depth 
requires a radiation transport analysis 

Analytical approximations are used in 
some codes to avoid time consuming 
Monte Carlo techniques 


Monte Carlo radiation transport 
techniques (e.g., ITS, EGS) provide 
better results 

- Empirical parameterization of 
Monte Carlo output provides 
fidelity of Monte Carlo radiation 
transport results with speed of 
analytical solution 

- Monte Carlo results are best for 
general solutions 


ENERGY DEPDSITBON (MpV g ' nm^i # electrons/source particle 


Radiation Transport 
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Internal Charging: Electrostatic Models 



* Generation-Recombination (GR) Model 


Microscopic model explicitly treats the radiation generated charge 
carrier pairs, field induced drift of carriers, and loss through 
recombination 

9E 

e— = P+~P--pj- (4) 
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* Radiation Induced Conductivity (RIC) Model 

Macroscopic model based on empirical relation ship between 
radiation dose and conductivity 


c IE 
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G FI model 


RIC model 


electron carrier generation 

deposition or RIC 


[Sessler et al. 2004] 
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Internal Charging: Electrostatic Models 



* Generation-Recombination (GR) Model 


Microscopic model explicitly treats the radiation generated charge 
carrier pairs, field induced drift of carriers, and loss through 
recombination 
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* Radiation Induced Conductivity (RIC) Model 

Macroscopic model based on empirical relation ship between 
radiation dose and conductivity 


dE _ 



dJp ■ d[(tr D +kP x )E] 
dx dx 



csda-range 





G FI model 


RIC model 


electron carrier generation 

deposition or RIC 


[Sessler et al. 2004] 
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DICTAT Internal Charging Code 


DICTAT model evaluates possibility of dielectric breakdown of insulating 
materials due to radiation charging 

- Single material parameterized by electrical conductivity (dark and radiation induced) 
dielectric constant, density, temperature 


J) Spacecraft charging: DICTAT Parameters - Mozilla Firefox 


Example: 1 mm Mylar exposed to FLUMIC (worst case) 
GEO electron environment 

• Resistivity 10 14 Q-m 

After 48.0 hours: Charging current= 7.0535E-13 Amps/cm 2 
E-max= 6.5170E+05 V/m Voltage= 464.5 volts 

The dielectric IS NOT liable to experience breakdown 
Maximum E lower than Breakdown field= 1.00E+07 V/m 

• Resistivity 10 16 Q-m 

After 48.0 hours: Charging current= 7.0535E-13 Amps/cm 2 
E-max= 1.4130E+07 V/m Voltage= 1.1678E+04 volts 

The dielectric IS liable to experience breakdown 
Maximum E higher than Breakdown field= 1.00E+07 V/m 

• Resistivity 10 18 Q-m 

After 48.0 hours: Charging current= 7.0535E-13 Amps/cm 2 
E-max= 3.5881E+07 V/m Voltage= 2.8170E+04 volts 

The dielectric IS liable to experience breakdown 
Maximum E higher than Breakdown field= 1.00E+07 V/m 


http://www.spenvis.oma.be/htbin/spenvis.exe/GTO 


SI 


1 Firefox prevented this site from opening a pop-up window, 


Options |j§ 


DICTAT Parameters 


Electron environment: | evaluate the FLUMIC model ata(B/Bo,L) location 


B/B 0 : 1.0 


L: 6.5 


Fraction of solar cycle: 1 0.0 


Fraction of year: 0.0 


Geometry and materials. Caution: the default values for 
the dielectric material parameters are given for 
reference only. The results of the simulation critically 
depend on the values of the input parameters. 


Geometry: 


Field of view 



Thickness [cm]: 
Temperature [K] : 

Density [g cm' 3 ]: 
Conductivity [Ohm 'm 1 ]: 
Dielectric constant: 
Breakdown el. field [V m" 1 ]: 
RIG dose rate factor k p : 
Delta : 

Activation energy [eV]: 


298.0 


1.0E-18 


3.0 


3.0E-17 




0.0 


Grounded at one surface 


[inner v| surface 


I aluminium 


Thickness [cm]:[£^_ 


Reset | [ Run | 


http://www.spenvis.oma.be/ 
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DICTAT Internal Charging Analysis 



• Simple first order evaluation of charging 
along trajectories through Earth's 
radiation belts using standard trapped 
radiation models, mitigation of charging 
by shielding 


• 10 mm thick plexiglass shielded by 

aluminum 0.0 cm to 0.2 cm thick 

- Peak AE-8 GEO radiation flux (constant) 


10 cm plexiglass (a=10 15 S/m) 
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Q IQ 20 3Q 40 SO 60 

Charging Time (hours) 


• Three 0.1 materials with conductivity 
varying from 10 16 to 10 12 S/m 

- AE-8 solar max electron flux 

- 250 km x 38,226 km x 0° inc, single orbit 

• Grey levels indicate relative threat levels 
for electrostatic discharge assuming 
materials typically suffer dielectric 
breakdown at electric fields in the range 
of 10 6 V/m to 10 7 V/m. 



o 


5 1 O 

Charging Time {hours) 


IS 
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DICTAT Internal Charging Analysis 


Electric field, surface voltage 
potential as function of time 

Allows quick screening for 
materials issues 

Simulation options limited by tool 

- Maximum time (48 hours) 

- Single material 

- Constant flux for user input 

- Radiation belt models for time 
variations in radiation flux 

- Only maximum E field is 
reported, no information on 
electric field as function of depth 





Electron Flux (#/cm2-5ec-sr-MeV) 



NUMIT ("numerical iteration") Codes 


• Originally developed by A.R. Frederickson et al (AFRL, JPL) 

• Tabata algorithm fits to Monte Carlo electron transport code for radiation dose, charge deposition 

• Multiple material conductivity, dielectric constant 

• Input electron spectrum, monoenergetic beams, fixed flux or time series of arbitrary duration (hours to years) 

• Output electric fields, charge density, currents, radiation dose rates as function of depth in material 
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3 MeV electrons on PMMA 



Charging beam 

• 3 MeV electrons 

• 0.01 nA/cm 2 

Material 

• PMMA (acrylic) 

• Z=6, A=12 

• ~10 cm thick 

• a ~ 1x10 16 S/m 

• k = 3.71 

• p = 1.19g/cm 3 


£ 

E 


100 

80 

60 
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E-Field Magnitude 




Charge Density 



110.0 

8.6 

7.2 

5.8 

4.4 
I 3.0 


SI 

J o 
1 o 

i\ 

o 


3.0 MeV 


max | +E| 


max |-E| 


4 6 8 

Time {hour) 


10 


12 


12 


v/m| 



3 MeV electrons on PMMA 



Charging beam 

• 3 MeV electrons 

• 0.01 nA/cm 2 

Material 

• PMMA (acrylic) 

• ~10 cm thick 

• Z=6, A=12 

• a ~ 1x10 16 S/m 

• k = 3.71 

• p = 1.19g/cm 3 



Time {hour) 
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3 MeV electrons on PMMA 



Charging beam 

• 3 MeV electrons 

• 0.01 nA/cm 2 

Material 

• PMMA (acrylic) 

• ~10 cm thick 

• Z=6, A=12 

• a ~ 1x10 16 S/m 

• k = 3.71 

• p = 1.19g/cm 3 


Beam off at 1 hour 
No effect since 
x~K8 0 /a=91 hours 
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3 MeV electrons on PMMA 



Charging beam 
• 3 MeV electrons 

? 

• 0.01 nA/cm 2 

,£ 

X 

Material 


• PMMA (acrylic) 

? 

,£ 
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X 

• Z=6, A=12 


• c ~ 1x10 14 S/m 

• k = 3.71 

• p = 1.19g/cm 3 
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#/cm2-sep-sr-M eV 


Geostationary Orbit 200 Day Time Series 


LANL/GEO: 2006lanl01 Smin.sav 
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Geostationary Orbit, 200 Days 



Charging current 

• Energy, flux from GEO ^ 
electron measurement! 

• LANL-01A SOPA +ESA X 


Material 

• PMMA 

• Z=6, A=12 

• ~2.5 cm thick 

• a ~ 1x10 17 S/m 

• k=2.00 

• p = 1.0 g/cm 3 
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Internal Charging Codes 


Model 

Type* 

Electron Transport 

Features 

Reference 

NUMIT (original) AFRLJPL 

A, MC 

Tabata, Monte Carlo 

l-D, single energy, a(E) models 

Frederickson et al., 1974, 
1977, 1980, 1983, 1993 

NUMITJPL 

A 

Tabata 

1-D, spectrum, 3 materials 

Jun et al. 2008 

NUMIT MSFC 

A 

Tabata 

1-D, spectrum, 5 materials 

Minow et al. 2007 
Jun et al 2008 

NUMIT Bethel University, 
AFRL 

A 

Tabata 

1-D, 1 energy, 3 materials 

Beeken and Mclver, 2010 

NUMIT SAIC 

A 

Tabata 

1-D, spectrum, radiation shield 

Davis et al., 2000, 2007 

ESA Deep Dielectric Charging 
Code (ESA-DDC) 

MC 

Range-energy 

relationship 

1-D, ct(E) models 

Soubeyran et al., 1993, 
1994 

DICTAT 

A 

Range-energy 

relationship 

1-D, radiation shield, ct(E) models 

Rodgers et al. 1999, 2000 , 
2003, S0rensen et al. 2000 

Moscow State University 

MC 

GEANT-3 

1-D 

Mileev and Novikov, 2004 

Xi'an Jiaotong University 

A 

1-D analytical solution 

1-D, ct(E) models 

Li et al. 2010 

Assessment Tool of Internal 
Charging for Satellites (ATICS) 

MC 

GEANT-4 

1-D panels over a 3-D spacecraft, 
a(E) models 

Zhongetal., 2007 

Multi-Utility Spacecraft 
Charging Analysis Tool 


1-D 

Flatta et al. 2009 
Cho et al. 2010 

JPL 

MC 

ITS 

3-D 

Katz and Kim, 2010 

Aerospace Corp 

MC 

GEANT-4 

3-D 

Lemon et al. 2010 


*A = analytical MC = Monte Carlo 
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Summary 



• Internal charging is a risk to spacecraft in energetic electron environments 

• DICTAT, NUMIT computational codes are the most widely used engineering 
tools for evaluating internal charging of insulator materials exposed to 
these environments 

• Engineering tools designed for rapid evaluation of ESD threats, but there is 
a need for more physics based models for investigating the science of 
materials interactions with energetic electron environments 

• Current tools are limited by the physics included in the models and ease of 
user implementation... .additional development work is needed to improve 
models 
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